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Abstract 

We report on the first non-perturbative calculation of the scalar Yukawa model in the single-nucleon sector 
up to four-body Fock sector truncation (one “scalar nucleon” and three “scalar pions”). The light-front 
Hamiltonian approach with a systematic non-perturbative renormalization is applied. We study the n-body 
norms and the electromagnetic form factor. We find that the one- and two-body contributions dominate up 
to coupling a « 1.7. As we approach the coupling a « 2.2, we discover that the four-body contribution rises 
rapidly and overtakes the two- and three-body contributions. By comparing with lower sector truncations, 
we show that the form factor converges with respect to the Fock sector expansion. 
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1. Introduction 

Solving quantum field theories in the non- 
perturbative regime is not only a theoretical chal¬ 
lenge but also essential to understand the struc¬ 
ture of hadrons from first principles. The light- 
front (LF) Hamiltonian quantum field theory ap¬ 
proach provides a natural framework to tackle this 
issue mm- A great advantage of this approach is 
that it provides direct access to the hadronic ob¬ 
servables. In the LF dynamics, the system is de¬ 
fined at a fixed LF time x'^ = t + z. The physical 
states are obtained by diagonalizing the LF Hamil¬ 
tonian operator. The vacuum in LF quantization is 
trivial. As a result, it is particularly convenient to 
expand the physical states in the Fock space. For 
example, a physical pion state can be written in 
terms of quarks (g), antiquarks (q) and gluons (g) 
as Itt) = \qq) + \qqg) -f \qqgg) -P ■ ■ ■. 

In order to do practical calculations, the Fock 
space has to be truncated. A natural choice, taking 
advantage of the LF dynamics, is the Fock sector 
truncation, also known as the light-front Tamm- 
Dancoff (LFTD) [5]. A number of non-perturbative 
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renormalization schemes have been developed based 
on the LFTD [3 [U 0 |^. Thus we arrive at a 
few-body problem and predictions can be system¬ 
atically improved by including more Fock sectors. 
The LFTD method is a non-perturbative approach 
in Minkowski space, which can be compared with 
other non-perturbative methods, e.g. Lattice quan¬ 
tum field theory in Euclidean space. Of course, 
this approach only works if the Fock sector expan¬ 
sion converges in the non-perturbative region. In 
practice, one can compare successive Fock sector 
truncations and check numerically whether the rel¬ 
evant physical observables converge. We will see 
that good convergence is achieved for the scalar 
Yukawa model in a non-perturbative regime with 
a four-body Fock sector truncation. Similar re¬ 
sults, though by a different method, were found in 
Refs. 1311] for the Wick-Cutkosky model [3]. 

We apply this approach to a scalar version of 
the Yukawa model that describes the pion-mediated 
nucleon-nucleon interaction. The Lagrangian den¬ 
sity of the model reads 

if = - m^l Ap -f ^d^ird^n - 

+ go\N\^7r + Sm^\N\I ( 1 ) 

where go is the bare coupling, Sm^ is the mass coun¬ 
terterm of the field N{x). It is convenient to intro- 
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duce a dimensionless coupling constant 


167rm^ 

For the sake of brevity, we refer to the funda¬ 
mental degrees-of-freedoms (d.o.f.’s) N{x) and 7r(x) 
as “scalar nucleon” and “scalar pion” held respec¬ 
tively. We also introduce a Pauli-Villars (PV) 
scalar pion (with mass fj-i) to regularize the ul¬ 
traviolet (UV) divergence [TO]. Then, a sector de¬ 
pendent method known as the Fock sector depen¬ 
dent renormalization (FSDR) developed in Ref. [S] 
is used to renormalize the theory. FSDR is a sys¬ 
tematic non-perturbative renormalization scheme 
based on the covariant light-front dynamics (CLFD, 
see Ref. m for a review) and Fock sector expan¬ 
sion. It has shown great promise in the application 
to the Yukawa model and QED [T2 US] . 

The scalar Yukawa model is known to exhibit 
a vacuum instability M- It can be stabilized 
by either adding the quartic terms 
and to the Lagrangian, or restricting the 

nucleon-antinucleon d.o.f. m The latter leads to 
the exclusion of the pion self-energy correction, 
sometimes referred to as the “quenched approxi¬ 
mation” . For the sake of simplicity, here we study 
this restricted version of the theory. Then the bare 
mass of the scalar pion becomes the physical mass, 
/io = 1^- It should be emphasized, though, that our 
formalism is capable of dealing with the (scalar) 
antinucleon d.o.f. The scalar nucleon and scalar 
pion d.o.f.’s generate non-perturbative dynamics at 
large coupling sufficient for our purposes. 

Previously, this model has been solved in the 
same approach up to three-body truncation (one 
scalar nucleon, two scalar pions) [3]. The results 
from the two- and three-body truncations agree at 
small couplings; yet they deviate in the large cou¬ 
pling region. Therefore, it is crucial to extend the 
non-perturbative calculation to higher Fock sectors. 
In this paper, we present the calculation of the four- 
body truncation (one scalar nucleon, three scalar 
pions). By comparing successive truncations, we 
can examine the convergence of the Fock sector ex¬ 
pansion. We presented a preliminary version of this 
work in Ref. |16) . 

We first introduce our formalism in the next sec¬ 
tion. The LF Hamiltonian field theory will be 
briefly mentioned and the non-perturbative renor¬ 
malization procedure will be explained. Then a set 
of coupled integral equations will be derived for the 
four-body truncation. In Sec. we present the nu¬ 


merical results, including the calculation of the elec¬ 
tromagnetic form factor. We conclude in Sec. 


2. Light-Front Hamiltonian Field Theory 

The LF Hamiltonian for the scalar Yukawa model 
is 


P- = 


■d±N+'m‘^\N\'^ + ^d±TT-d±TT 


+ - 5o|IVp7r - 5m^\N[" 


. rE+=0 


■ ( 2 ) 


The physical states can be obtained by solving the 
time-independent Schrodinger equation 


p- Ip) = Ip), 

p+ 


(3) 


where pj_ and p'*" are the transverse and longitudi¬ 
nal momentum, respectively. Thanks to boost in¬ 
variance in the LF dynamics, we can take pj_ = 0 
without loss of generality. 

The system is solved in the single-nucleon sector. 
The state vector admits a Fock space expansion. 


\p) = Y. 

n 


Pji '^n (I^1_L 5 5 ‘ ’ ' ^n_L i ) 


X |fci_L,a;i,---fe„_L,a:„), (4) 


where Xi = 


p+ 5 


and 


Dn = 2(27r)^d(^) (feij_ H- kni_)5[xi H-- 1) 

{2TTY2xi 


The n-body Fock state |fei_L, cci, • • • ,kn±,Xn) con¬ 
sists (n — I) scalar pions and 1 scalar nucleon. We 
use the last pair {kn±,Xn) to denote the momen¬ 
tum of the scalar nucleon, pn, known as the LF 
wave function (LFWF), is a boost invariant. The 
LFWFs are normalized to unity, J2n -In = !> where 


= 


(n-l)! 


-D. 


^n(I^l_Lj^l; ‘ ‘ ' kjij_ ^ Xrfi^p ) 


( 5 ) 

is the probability that the system appears in the 
n-body Fock sector. In the scalar Yukawa model, 
these quantities are regulator independent, in con¬ 
trast to more realistic theories such as Yukawa and 
QED. Note that pi = ^/li is a constant. 
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Figure 1: The diagrammatic representation of the system of equations in the four-body truncation. 



^2 Ti r-i Ti Ti 


Figure 2: The perturbative expansion of the two-body vertex function. The solid lines represent the scalar nucleons; the 
dashed lines represent the scalar pions; the double lines represent the dressed nucleons. 


It is convenient to introduce the n-body vertex 
functions, 


1 1_L : ^n— 1: P ) — 

(-Si,... ,n—1 P )t/’n (^1_L ) it'1; ‘ ‘ ‘ ^ ^n_L : ? P ) (6) 

for n > 1 and Ti = (rn? — p'^)ipi, where 



i—ii 


is the invariant mass squared of the Fock state, and 
Pji Ui — Oj 1) mass of the i-th scalar pion. We 
have suppressed fc„_L and a;„ in r„, by virtue of the 
momentum conservations ki± -|- k 2 ± + ■ • • fc„_L = 0, 
xi + X 2 + ■ ■ ■ + Xn = ^- For simplicity we will also 
omit the dependence on in r„ for the ground 
state p^ = 

Written in terms of the vertex functions F, 
Eq. (§ can be represented diagrammatically using 
the LF graphical rules [T71 [TS] (see Ref. [TT] for a 
review). Figure]^ shows the diagrams for the four- 
body truncation. 

The two-body vertex function r 2 plays a partic¬ 
ular role in renormalization. It comprises all radia¬ 


tive corrections allowed by the Fock sector trunca¬ 
tion, including the amputated vertex V 2 ,{ki,k 2 ,p) 
and the self-energy E((p — fci)^) (see Fig. [^: 

F2(fci_L,Xi;p^) = Z{{p- kif )V3{ki,k2,p)\/h, 

(7) 

where the function Z{q^) = (l — ^ is a 

generalization of the field strength renormalization 
constant Z = {l - = h. Note the 

presence of the scalar pion spectator, which means 
that in the n-body truncation, the self-energy cor¬ 
rection in the expression for F 2 is the (n — l)-body 
self-energy. 

The dependence of renormalization constants on 
the Fock sector is a general feature of the Fock 
sector expansion. We use gon and Srn^ to de¬ 
note the bare coupling and the mass counterterm 
from the n-body truncation, respectively. Accord¬ 
ing to the LSZ reduction formula, the physical cou¬ 
pling g = Tfi = '/ZV^{ki,k 2 ,p)^/h. Here 
means that V 3 is evaluated at the renormalization 
point, the physical mass shell si = mZ => k\j_ — 
— {l — xi)pL 3 — x\mZ = k^ 3 _. These relations provide 
the on-shell renormalization condition [susma, 

( 8 ) 
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Figure 3: The self-energy correction, loop correction E plus 
mass counterterm Sm?, expressed in terms of the two-body 
vertex function r 2 . Note the external lines are amputated. 


Here the Fock sector dependence is shown explic- 
itly. For example, F^ ^ represents the two-body ver¬ 
tex function found in the n-body truncation. Note 
that k*'j_ is negative, which means Eq. (|^ has to 
be imposed through analytic continuation. 

The two-body vertex function F2 also provides a 
non-perturbative means to calculate the self-energy 
correction (see Fig. [^. Following the LF graphical 
rules, the self-energy in the n-body truncation is. 


'r(")W2 f dxi gon 

. W J (27r)3 Jo 2 x 1 ( 1 -Xi) 


r^^'^"(fci_L,a;i,fc 2 _L,X 2 ) = 




503(61) r2"(fei-L:a^i) 


(1 - Xi)(si - m?) 


'5o2 


1 

E 

J3=0 


(- 1 ) 


,J 3 


1—XI—X2 

d^6_L f 


dX 3 


(27r)3 J 2x3(1 - Xi - X 3 )(l - Xi - X 2 - X 3 ) 
0 


1 F^^-^yfci_L,Xi,fc 3 _L,X 3 ) 
S123 - (si3 - 


+ (1 O 2) 

( 11 ) 


* _ 


l +6 


where 61 = X2/(1-Xi), s^2 = xi ■ X2 • 

and Z(2) comes from combining the two-body self¬ 
energy corrections. We have included the PV 
scalar pions (j = 1) in the equations along with 
the “physical” pions {j = 0 ). As mentioned, 
502, ^m2,5o3, are sector dependent renormal¬ 
ization “constants” obtained from the two- and 
three-body truncations [ 5 ]. In fact, (703 depends 
on the momentum fraction x, which is a manifes¬ 
tation of the violation of the Lorentz symmetry by 
the Fock sector truncation [ 12 ]. F3 is an auxiliary 
function that satishes the integral equation. 


Note that in our formalism the state vector in 
Eq. rather than its one-body component, is nor¬ 
malized to unity. So according to the definition of 
the self-energy, the one-body LFWE ipi = 1/J1 is 
excluded from F2 in the above expression. Then 
the mass renormalization condition in the on-shell 
scheme implies Srri^ = S*^")(m^). 


r3'’Mfel_L:a;i,fc2_L,X2) = 




503(62)r^6fc2.L,X2) 
(1 - X 2 )(S 2 - m 2 ) 


1 


+ 50^2 

33^0 


X 


I — X1—X2 

d^fcs-L f _dX3_ 

( 27 r )3 J 2x3(1 - Xi - X 3 )(l - Xi - X 2 - X 3 ) 
0 


As mentioned, the system of equations for r2-4 
resulted from truncating Eq. ^ to at most four- 
body (one scalar nucleon and three scalar pions) 
are shown in Fig. After substituting F4 into the 
second equation and applying the renormalization 
condition Eq. ([^, the system of equations becomes 


r3'^"(^lJ_>®l’^3_L,X3) 


’123 \ '’13 

^ r^"-^yfc2_L,X2,fc3_L,X3) 

S23 - m2 


( 12 ) 


where 62 = Xi/(1 - X2), q*^ = m^ - (1 - Xi - 


'' (1 - Xi)(si - m2) 


E(-i) 

12 = 0 


r^^'^"(fcl-L,Xi,A;2_L,X2) 


Si 2 — m^ 
r3'’6fel_L>2:i,fc2_L,X2) 


’12 


(10) 


X 2 )(s^ 2 -m 2 ), ^23 = + 


kix+t^r 


(klx+k 2 ±+k 3 ±_)'^ +r, 

I — X1—X2 — X3 


Note that Eq. ( 10 1 can be 


X3 


f <i‘^k 2 ± 

f dx 2 503(61) 

Eq. (12 

J ( 2^)3 j 

2 x 2(1 - Xi - X 2 ) 



3. Numerical Results 

We employ an iterative procedure to solve 


Eqs. (10 - 12 ). The momenta are discretized on cho¬ 


sen grids in the transverse radial and angular co¬ 
ordinates as well as in the longitudinal coordinate. 
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Figure 4: The Fock sector norms /i _4 as a function of the 
PV mass /ri for a = 1.0 (top) a = 2.0 (bottom). Results 
evaluated on different grids are shown. 


where the grid sizes are controlled by the number of 
abscissas, fV^ad, -^ang, and Nn^. Then the integrals 
are approximated by the Gauss-Legendre quadra¬ 
ture. We start with an initial guess of the vertex 
functions and update them iteratively, until reach¬ 
ing a pointwise absolute tolerance max{|Ar|} < 
10“'^. We solved the system at to = 0.94 GeV, /to = 
0.14 GeV. The numerical results are obtained using 
Cray XE 6 Hopper at NERSC. 

Figure]^ plots the Fock sector normalization fac¬ 
tors In (see Eq. §) as a function of the PV mass 
/ii for two selected coupling constants. It shows 
that for sufficiently large grids, /„ converge as /ti 
increases. However, for a fixed grid, increasing /ti 
would increase the numerical error while decreasing 
the systematic error introduced by the finite regula¬ 
tor, as larger /^i requires more coverage in the UV 
hence larger grid size. A PV mass = 15 GeV 
suffices for our purposes here. 

There exist two critical couplings at ac ~ 2.6 and 
a(, « 2.2. In the two-body truncation, one finds the 



Figure 5: The Fock sector norms /i —4 as a function of 

the coupling constant a. Results are evaluated on the grid 
TVifx = 41, A^rad = -^aiig = 20, with a PV mass pi = 15 GeV. 


bare coupling. 


1 1 
9^ 5o2 




. \m / 

\m /. 


where/(A) = dxx{l — x)/{{l — x)X‘^ + x^) and 
/(A —>■ oo) = 0. If the physical coupling constant 
a > ac = 7r//(/io/TO), the two-body bare coupling 
<702 diverges at some finite PV mass. Such a sin¬ 
gularity in (/02 (known as the “Landau pole” in a 
similar case in QED) propagates from the two-body 
truncation to the four-body truncation via <702 used 
in the FSDR. At a = «(,, the determinant of the 
Hamiltonian in the three-body truncation crosses 
zero. Similarly, this singularity propagates from 
the three-body truncation to the four-body trun¬ 
cation and the iterative procedure in the four-body 
truncation diverges at a > a(,. 

Figure shows the contribution of each Fock 
sector in the four-body truncation for couplings 
up to a = 2.12. A natural Fock sector hierar¬ 
chy Ii > I 2 > I 3 > I 4 can be observed, up to 
a « 1.7. Beyond a « 1.7, J 4 exceeds I 3 and be¬ 
gins a steep climb with increasing a. Meanwhile, I 2 
turns over and starts to fall. The net effect is that 
I 4 exceeds I 2 and I 3 at about a « 2.1. Clearly, 
as we approach «(,, dramatic changes in the /„’s 
are emerging and it appears that the Fock space 
expansion breaks down. Nevertheless, the lowest 
sectors \N) + \NTr) are observed to dominate the 
Fock space up to a « 2.0, where these two sectors 
constitute 80% of the full norm. 

Figure compares the Fock sector norms from 
the four-body truncation with their counterparts 
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Figure 6: Comparison of the Fock sector norms Ii (top) 
and I 2 (bottom) from successive two-, three- and four-body 
truncations. 


from the two- and three-body truncations. The re¬ 
sult suggests a convergence as the number of con¬ 
stituent bosons increases, especially for the cou¬ 
pling below a ~ 1.0. Note that the one-body norm 
Ii changes little from the three-body truncation to 
the four-body truncation, even around a ~ 1.7. 

The obtained LFWFs are now available for com¬ 
puting physical observables. Here we consider the 
elastic electromagnetic form factor for photon cou¬ 
pling to the scalar nucleon, which is obtained from 
the matrix element of the “-I-” component of the 
current (see Fig. [^, 

(p + glJ+(0)lp} = 2p+F(Q^), (13) 

where g~^ = 0, = —g^ = > 0. In LF dynam¬ 

ics, the form factor obtains the form [19]: 

n 

X 'tpniki±,Xi ) (14) 




(a) 


(b) 



Figure 7: The various Fock sector contributions to the elec¬ 
tromagnetic form factor for the scalar nucleon (solid lines). 
The dashed lines represent the internal scalar pions. The 
wavy lines represent the external photons. 


where fc'j_ = ki± — Xiqj_, {i = 1, 2, • • • , n — 1), for 
the spectators and = fc„_L + (1 — a;„)qj_ for the 
struck parton. Figure shows the form factor for 



Figure 8: The elastic electromagnetic form factor F(Q^) in 
the four-body truncation for couplings a = 0.2, 0.4, 0.8, 
1.6, and 2.1. The numerical results (symbols) are fitted by 
Eq. ( |15| | (lines). 


some selected couplings. In the limit of —)■ 0, 

F(0) = 1, consistent with the charge conservation; 
in the limit of —)■ oo, F{oo) = Ii, represent¬ 

ing a point-like charge. The form factors can be 
approximated by 


Fm 


1 + c I\Q^ 

l + cQ2 


(15) 


Figure [^compares the form factors obtained from 
the two-, three- and four-body truncations for two 
selected couplings. The three- and four-body trun¬ 
cation results show good agreement even at the 


6 































- PV mass 15GeV 

n pi I I I I I I I I I I I I I I I I I I I I 

■ 0 25 50 75 100 

Cf [GeV^] 

Figure 9: Comparison of the form factors calculated in the 
two-, three and four-body truncations at a = 1.0 (top) and 
a. = 2.0 (bottom)^ The three- and four-body form factors 
are fitted by Eq. 

non-perturbative couplings, suggesting a reasonable 
convergence with respect to the Fock sector expan¬ 
sion. 


4. Discussion and Conclusions 

We solve the single-nucleon sector of the scalar 
Yukawa model in light-front dynamics within a 
four-body (up to one scalar nucleon and three scalar 
pions) Fock sector truncation. Fock sector depen¬ 
dent renormalization is implemented. The coupled 
system of linear integral equations is derived and 
solved numerically. The numerical study of the 
Fock sector norms suggests that up to a « 1.7 
the system is dominated by the lowest Fock sec¬ 
tors. By comparing the form factors from succes¬ 
sive Fock sector truncations (two-, three- and four- 
body) , we find that the Fock space expansion of the 
form factor for the scalar nucleon converges as the 


number of scalar pions increases even in the non- 
perturbative region. 

Solving the one-nucleon sector is also the first 
step for the study of the two-nucleon sector - a 
bound-state problem, which has been extensively 
studied in various approaches (see, e.g., [5D] and 
the references therein). However not all these ap¬ 
proaches are from first principles. In our approach, 
the two-nucleon sector obeys similar integral equa¬ 
tions. The bare couplings and the mass countert¬ 
erms, according to FSDR, are already provided by 
the one-nucleon sector (up to three dressing pions). 
Therefore, our approach allows a systematic study 
of the theory with a non-perturbative renormaliza¬ 
tion. 

This calculation demonstrates that the light- 
front Tamm-Dancoff, equipped with the Fock sec¬ 
tor dependent renormalization, is a general ab initio 
non-perturbative approach to quantum field theo¬ 
ries. While the solution of the scalar Yukawa model 
may be useful for, e.g., chiral effective field theory 
studies, this approach has also been applied to more 
realistic field theories, including the Yukawa model 
(truncation up to one spinor and two scalars) [13] 
and QED (truncation up to one electron and two 
photons) [ 3 . In these theories, the vertex func¬ 
tions also diverge, in contrast to the scalar Yukawa 
model. However, after the renormalization, the 
physical observables converge as expected. Never¬ 
theless, the study of the higher Fock sector expan¬ 
sion in these models is in principle similar to the 
current one, which indicates the potential of this 
approach as an alternative to other first-principle 
methods, e.g. the lattice gauge theory, especially 
in the study of hadron structures. 
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